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ABSTRACT 

In this paper we present a time-varying sliding mode control (TVSMC) technique for 
reusable launch vehicle (RLV) attitude control in ascent and entry flight phases. In ascent 
flight the guidance commands Euler roll, pitch and yaw angles, and in entry flight it 
commands the aerodynamic angles of bank, attack and sideslip. The controller employs a 
body rate inner loop and the attitude outer loop, which are separated in time-scale by the 
singular perturbation principle. The novelty of the TVSMC is that both the sliding surface 
and the boundary layer dynamics can be varied in real time using the PD-eigenvalue 
assignment technique. This salient feature is used to cope with control command saturation 
and integrator windup in the presence of severe disturbance or control effector failure, which 
enhances the robustness and fault tolerance of the controller. The TV-SMC ascent and 
descent designs are currently being tested with high fidelity, 6-DOF dispersion simulations. 
The test results will be presented in the final version of this paper. 
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1. Introduction 


Flight control of both current and future reusable launch vehicles (RLV) in ascent and 
descent modes involves attitude maneuvering through a wide range of flight conditions, wind 
disturbances, and plant uncertainties including aerodynamic surfaces and engine failures. The 
baseline RLV flight control system that was designed for the X-33 technology demonstration 
sub-orbital launch vehicle employs a variable structure PID control law [1] with gain schedulmg. 
This requires four gains per channel that are looked up from a table as a function of relative 
velocity during ascent flight. Depending on the flight trajectory, each gain table can have as many 
as 25 values, so potentially 300 gain values must be stored in the on board computer for nominal 
flight. In case of one engine failure, or Power Pack Out (PPO), alternate sets of gain tables are 
used, depending on the flight time when the failure occurred. Provisions are made for 25 possible 
PPO times, or 25 sets of PID tables. This amounts to another 7500 values to be stored, or a total 
of 7800 values to provide gains for the nominal and engine failed cases during ascent flight. The 
reason for so many gain tables is because the control system design relies on linear analysis and 
perturbation theory at specific design points along the trajectory. This method is well established 
and has been used in many launch vehicle control system designs. Satisfactory control 
performance is ensured as long as the vehicle performance and operating conditions are relatively 
close to the design points. A robust flight control algorithm that would accommodate different 
trajectories and aerodynamic surface and engine failures without gain scheduling would be an 
improvement over the RLV current flight control technology. Sliding Mode Controller (SMC) is 
an attractive robust control algorithm for the RLV ascent and descent flight controller designs 
because of its inherent insensitivity and robustness to plant uncertainties and external 
disturbances [2-4]. Such a robust controller [5,6] would reduce risk and drastically decrease the 
amount of time spent in pre-flight analysis, thus reducing cost 

In this work we present a time-varying sliding mode controller (TV-SMC) design technique. 
The RLV fixed-gain (FG) SMC two-loop structure (Fig. 1), which is presented in the work [5,6], 
is employed in this work for the TV-SMC design. 
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Fig. 1 The RLV SMC two-loop structure 

In the outer loop, the kinematics equation of angular motion is used with the outer loop 
SMC to generate the angular rate profiles as virtual control inputs to the inner loop. In the inner 
loop, a suitable inner loop SMC is designed so that the commanded angular rate profiles are 
tracked. The inner loop SMC produces roll, pitch and yaw torque commands, which are allocated 
into end-effector deflection commands. Multiple time scaling (multiple-scale) is defined as the 
time-constant separation between the two loops. That is, the inner loop compensated dynamics is 
designed to be faster than the outer loop dynamics. The resulting multiple-scale two-loop SMC, 
with optimal torque allocation, causes the angular rate and the Euler angle (ascent) or 
aerodynamic angle (descent) tracking errors to be constrained to linear de-coupled homogeneous 
vector valued differential equations with desired eigenvalues. 

The TV-SMC design technique is developed for improved robustness in the presence of 
control effector saturation and tolerance of control effector failure. Any partial or complete 
failure of control actuators and effectors can be inferred from saturation of one or more 
commanded control signals generated by the controller. The saturation causes a reduction in the 
effective gain, or bandwidth of the feedback loop. The controller cannot tell nor does it care 
whether control command saturation is due to wind gust or a stuck effector, because it only 
computes the required torque to eliminate the tracking error. A truly adaptive controller should 
be able to do this even if the vehicle has lost significant control authority due to control effector 
failure. It is up to the combination of the control allocation and the control effectors to realize the 
required torque. In order to maintain stability, the bandwidth of the nominal (reduced-order) 
system will be reduced accordingly using a time-varying bandwidth PD-eigenstructure 





assignment technique [7,8]. The presented TV-SMC fault-tolerant technique automatically 
handles momentary saturations and integrator windup caused by excessive disturbances, 
guidance command or dispersions under normal vehicle conditions. 

The TV-SMC ascent and descent designs are currently being tested with high fidelity, 6- 
DOF dispersion simulations. The test results will be presented in the final version of this paper. 


2. RLV Attitude Control Problem Formulation 

The rigid body equations of motion for an RLV is given by (l)-(5) 

(j 0 +Aj)d> = -£2(J 0 +AJ)y+(l+E)r+d, (i) 
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The control problem for the RLV in ascent and descent modes is to determine the control torque 
command vector T such that the commanded orientation angle profiles y c are robustly 
asymptotically followed in the presence of bounded disturbance torque d, the RLV inertia 
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variations AJ and aerodynamic surface and engine failures that are described by the uncertain 
matrix AD(.) , E = AD(.)D* , D*(.) = D(.) T [d(.)D(.) t Y . D(.)is a nominal sensitivity matrix. 

3. Summary of the Smooth Multiple-loop Fixed-Gain Sliding Mode Control 

The smooth multiple-loop fixed-gain sliding mode controller is designed as follows [5,6]: 

• The outer loop smooth fixed-gain SMC generates body rate commands 

<y c =R"'(r)[f c +K l rJ+R' I (r) K 2<7 

(T = y e +K l \y e dT, creR\ y t = Y e - Y V) 

0 

or 

t 

co c =R- l (y) ^+(K 1 +K 2 )r e +K 2 K l f y e dt (8) 

0 

that provide the following outer loop compensated dynamics (given the 0) c profile is 
tracked perfectly in the inner loop): 

f'+iKt+K^K+iK^y^Q. ( 9 ) 

• The inner loop smooth fixed-gain SMC generates control torque commands 

T = J 0 d) c + J 0 K 3 6i e +nJ 0 <y+J 0 K 4 s (10) 

/ 

s = (o e +K 3 \o) e dr, seR 3 , 0 ) e = 0 ) c - 0 ) (11) 

0 

or 

/ 

T = 3 0 d) c + Q3 0 o) c + (J 0 K 3 -Q3 0 + J 0 K 4 K +J 0 K 4 K 3 f 0) t dr (12) 

0 

that provide the following inner loop compensated dynamics (given the T profile is 
allocated perfectly into commands to actuator deflections): 

d> e + (K 3 +K 4 )f r + (K 4 K 3 )n=0 (13) 
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• Writing inner and outer loop tracking error equations in a damping factor/natural 
frequency format 

x+2ga„x+£i>Zx = 0 ( 14 ) 

it is easy to calculate elements of diagonal matrices K, K 2 , K 3 and K 4 to provide given 

damping factors (usually all are equal to 1.1 or so) and natural frequencies that provide 
for a sufficient time-scale separation between the control loops to eqs. (9) and (13) for 
inner and outer loop compensated tracking dynamics. 

• Commands to the actuator deflections can be calculated as follows 


|A = d(.)ti 


t - +m 0 o) c + ( j 0 k 3 - j 0 k 4 X +/»2 j«k 4 k , f 


(15) 


where S c0 is a bias command that should trim the RLV, D(.) is a sensitivity matrix (for a 


feedforward control allocation case). 


4. Smooth Multiple-loop Time-Varying Sliding Mode Control 

Motivations. The elements of the command torque vector (12) are limited by physical 
abilities of the RLV flight control system that implies the following limits to elements of the 

vector T : 


a,QZT t <*, (•), i = L3. 

(16) 

a\Q^T^b]Q, * = U 
at 

(17) 


where a,Q. b,Q, a,Q, b]Q depend on the RLV current flying conditions including Mach 
number, dynamic pressure and trim conditions. Designing the time-varying SMC inequalities 
(17) could be originally out of consideration. 

The components of the vector T that depend on body rate tracking profile, eo c , and its 
derivative, d) c , usually have the largest amplitudes among the others. So, the following idea is 
proposed: if the inequalities (16) and/or (17) are about to be violated we have to start reducing 
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magnitudes of the terms co c and d) c to prevent actuators from saturation. It will be achieved via 

real time adjusting bandwidths of inner and outer flight control loops. This procedure requires 
application of time-varying linear control technique [7,8] that is incorporated into the time- 
varying sliding mode controller design. 

The following outer, (O c , and inner, T , continuous, time-varying sliding mode control 
(TV-SMC) laws are designed 

a) e = R" 1 (y)[? e +K 1 (0rJ+R' 1 (r) K 2(0o- < 18) 

T = J 0 <2J C ■*" ^ (O^e (^) 

where the outer, < 7 , and inner, s , sliding quantities are identified 

O’ = n + (K , (0 - K 2 (r)" 1 K, (r))f Y.dz , <reR 3 , Y,=Y C ~Y ( 20 ) 

0 

t 

s = C 0 e + (K i (t)-K 4 (tr l K i W)\G) e dT, se R 3 , Q) e =( 0 c - 0 ) (21) 

0 

Remark. One can observe an apparent difference between the eqs. (20), (21) and eqs. (7), (1 1). 
These equations become identical if K ; = const and K 3 = const . 

The diagonal elements of the matrices K 2 (t) and K 4 (t) determines the time-varying 
thickness of the boundary layers of the outer and inner loops, respectively. Once K 2 (t) and 
K 4 (t) are fixed, K x (t) and K 2 (t) define the dynamics of the sliding modes of the outer and 

inner loops, respectively. When these controller gains are chosen in a certain manner as defined 
below, they can be used to define the desired time-vaiying bandwidths of the outer and inner 
loops. These bandwidths can be adjusted in real time to prevent actuators from saturation. The 
corresponding bandwidths or natural frequency of the compensated loop dynamics are to be 
decreased if the actuator command is close to saturation and are to be increased otherwise. 
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5. Controller Design and Stability Assessment: Nominal Case 

In the nominal case where E = 0, AJ = 0, d = 0 , the mathematical model of the RLV in 
eqs. (l)-(5) is completely known. For the outer loop, substituting eqs. (18) and (20) into eq. (2) 
and differentiating eq. (2) we obtained the compensated outer loop error dynamics 

t 

^, + (K 1 (0+K 2 (0)? , . + ^i(t)Ka(0-K 1 (/)^ = R at. Y. = \Y MW ( 22 > 

o 


Assume that the transient response in the inner, "faster", loop is stable and (Q-(Q C in eqs. (2) and 
(18). Then the tracking error S) t is stabilized in the inner loop at zero. Assuming <B t = 0 the 
equation (22) can be rewritten in a homogeneous format 

f. + (K t (0 + K 2 (0)f t + (Ki(/)K 2 (r)-Ki(/)fc = 0 (23) 


Since the matrices K x (t) and K 2 (t) are diagonal, eq. (23) can be rewritten in a scalar format 

h + (K, (0 + *» (0)fw + {K (0 (0 - ki (OK = 0 , Vi = u (24) 

The series D-eigenvalues of LTV differential equations (24) can be identified as follows 

K (0 = ~ k u ( 0 . (0 = - k 2i( { )> Vi = 1,3 . (25) 


In order to assign values for the bandwidth (cut-off frequency or natural frequency) and damping 
factor for the compensated dynamics of the outer loop eqs. (24) are rewritten in a format 


Yei + 
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+ a>l (0 Y„- =0, >1 Vi = 1,3. 
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Then the so-called "PD-eigenvalue" can be identified as follows: 

A,(o=(-6+ 1 /iM)>a') < 27 > 
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Note that for time-vaiying co^t ) , the relations between the sign of the radicals and the index of 
A, must be maintained. Then the SD-eigenvalues are related to the PD-eigenvalues as follows 

4 (?) = K (0 =(-£+ V6 2 (0 = ~ k u (0 > V/= 1>3 (29) 


or 




(30) 


and 


= = MO (31) 


or 


^(0 = (^+^ 2 - 1 )® n i ( 0 - 


Ml 

® ni (0 


(32) 


Thus, given £ and <y ni (?) the corresponding time-vaiying coefficients in the sliding surface (8) 
can be computed using eqs. (30) and (32). 


Remarks. 

1. The outer loop dynamics are obviously decoupled and the stability is guranteed by the PD- 
eigenvalues for £ >1 and <0„(O > 0, V? >? 0> which can be achieved by the choice of 

diagonal time-varying matrices K { (?), K 2 (?) according to (30) and (32). 

2. The gain matrices K l (t),K 2 (t) must be selected in accordance with (30) and (32) with 
£ > i in order to to make a - dynamics faster then y t -dynamics. The time scale separation 
increases with £.Note that (30) and (32) are for constant £ only, and as such £ should 
not be adjusted in real time. The PD-eigenvalue synthesis formulas for time-varying £(?) 
is available, but it is much more complex. 
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The inner loop compensated error dynamics are described by the eqs. (1), (3), (19) and (21), 
which can be combined as 

/ 

a. +(k j (o+M>M +(k,( 0*.(')-*.('))3. =o. ®. = f «>.<*)* < 33 > 

0 


Eq. (33) can be rewritten in a scalar "damping factor-natural frequency" format 


co el + 


2|a CO- 


M2 

£ ni (0 


(O c > + (0 co ei — 0 , > 1 Vi — 1,3 . 


(34) 


Eqs. (33) and (34) is a full anology to eqs. (24) and (26). So, the following coefficients can be 
identified by analogy 


*„«)=(£ < 3s ) 

*. w -(A + ^r^ i( (36> 

Now, given f, and d) ai ( t ) the corresponding time-varying coefficients in the sliding surface (21) 
can be computed concerning eqs. (18) and (19) Vi = 1,3 . 


Remark. The same comments for K t (0,K 2 (0 following (30) and (32) are applicable to 
K 3 (t), K A (t ) in relation to (35) and (36), and the s - dynamics and -dynamics. 


6. Stability Assessment: Perturbed Case 

In the perturbed case assume E* 0, AJ*0, d± 0 , unknown but bounded, i.e. the math 
model of the RLV in eqs. (l)-(5) is not completely known. For the outer loop, again assume that 
the transient in the inner "faster" loop is stable and 0) = 0) c in eqs. (2) and (18). Then the 
perturbed case coinsides with the nominal case, since eqs. (2), (4) and (5) do not contain 
uncertaintes and disturbances. 
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The inner loop error compensated dynamics are described by the eqs. (1), (3), (19) and (21). 
In the perturbed case this is 

a, + (K, (o + k< (t))8} e + ( k, (t)K < (o - k> (oK = m A A a) ov 

where F(a) e ,d) e ,ti> e ,t) depends on perturbation terms. The closed-loop stability of the perturbed 
system can be justified by the well-known results from the Lyapunov analysis [9]. 

Theorem. Given a perturbed linear time-varying dynamic system 

x = A(t)x+F(xj) ( 38 ) 

where F{x,t ) < 5 is a bounded, nonvanishing perturbation, i.e. F(0,t ) * 0, and the matrix A(t) 
is Hurwitz, i.e. there exists a pair of positive definite matrices (P(t),Q(t)) satisfying the 


Lyapunov equation 


P{t ) + P{t)A{t) + A T ( t)P(t ) = ~Q{t) 

(39) 

such that 


Cl |lx(0|r £ x T (t)P(t)x(t) < c 2 |x(r)| 2 

(40) 

^(oeC'MO^-^IKof 

(41) 

p(0^(0l^c 4 |x(r)fl 

(42) 

for all t > t 0 , where flv| = -Jv T v . For any 0 < 0 < 1 , define 



(43) 

If F(x, t)<.5 <S a , then every trajectory of the perturbed system is exponentially, ultimately 
bounded in the sense that for any r > 0 . and ||x(r 0 )| £r/k, there exists a t { >t 0 such that 

lx(r,r 0 )l<%(t 0 )||e- A< '- ) , t 0 <,t<t x . 

(44) 
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|*(Mo)| 1 ^*i» 


(45) 


where 


k = 



( 1-^3 
2 c 0 



(46) 


7. Summary and Conclusions 

The time-varying sliding mode controller (TV-SMC) design algorithm is developed for the 
2 nd generation reusable launch vehicle in ascend and descend modes. In order to maintain 
stability, the bandwidth of the nominal (reduced-order) system is reduced accordingly using a 
time-varying bandwidth PD-eigenstructure assignment technique. The presented TV-SMC fault- 
tolerant technique automatically handles momentary saturations and integrator windup caused by 
excessive disturbances, guidance command or dispersions under normal vehicle conditions. The 
TV-SMC algorithm has been coded for the X-33 technology demonstration vehicle in ascent and 
descent modes. The results of simulations of the designed TV-SMC will be presented and 
discussed in the final version of the paper. 
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